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Diluted Josephson-junction arrays in a magnetic field: phase coherence and vortex 

glass thresholds. 
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The effects of random dilution of junctions on a two-dimensional Josephson-junction array in a 
magnetic field are considered. For rational values of the average flux quantum per plaquette /, the 
superconducting transition temperature vanishes, for increasing dilution, at a critical value xs(f), 
while the vortex ordering remains stable up to xvl > xs, much below the value x p corresponding to 
the geometric percolation threshold. For xvl < x < x p , the array behaves as a zero-temperature 
vortex-glass. Numerical results for / = 1/2 from defect energy calculations are presented which are 
consistent with this scenario. 



dimensions, vortex 
perconducting film: 



74.50+r, 74.60.Ge, 64.60.Cn 



Vortex glass states in disordered three-dimensional 
superconductors have been the focus of much recent 
interestEnj. In the absence of screening, they are believed 
to have a true superconducting phase, with vanishing lin- 
ear resistivity, at finite temperatures. By contrast, in two 
ass modclsa and experiments on su- 
show that vortex glass order is de- 
stroyed at any finite temperature with a nonzero but ex- 
ponentially small resistivity. This zero-temperature vor- 
tex glass can be characterized by a thermal correlation 
length exponent ut which determines, for example, the 
current density scale, J n i ~ r r 1 +^r ^ wnere nonlinear be- 
havior shows up in the current-voltage characteristicsB'El. 
Recent_cstimates give Vt ~ 2 for various vortex glass 
modelsa. 

Randomly diluted Josephson junction arrays (JJ 
have been used to model disordered superconductors!! 
In zero field, the superconducting transition temperature 
vanishes at the percolation thresholds x p , where x is the 
concentration of diluted junctions. For x > x p there 
are only uncoupled finite clusters and long-range phase 
coherence is destroyed. At x p , the infinite percolating 
cluster shows, up in the scaling behavior of the dynamic 
conductivityQ and nonlinear resistivityO. In the presence 
of an external field, a diluted JJA is an experimentally 
controllable model to investigate phase coherence and 
vortex glasses in two dimensions. For rational values of 
the flux quantum per unit cell /, an ordered (x = 0) JJA 
has a ground state consisting of a periodic pinned vor- 
tex lattice, with additional discrete symmetries resulting 
from commensurability effectsuiJ. The melting of this 
vortex lattice at a temperature Tvl, driven by domain- 
wall excitations, competes with the superconducting 
transition at T$ driven by the Kosterlitz-Thouless vortex 
unbinding. For / = 1/2, these transitions either coincide 



or have very close transition temperatures^, Tvl T$. 
Similar behavior is expected for other low rational values 
of /. In presence of random-dilution disorder, two natu- 
ral questions arise: (i) are there two dilution thresholds, 
xs and Xvl, for phase coherence and vortex lattice or- 
der respectively? Does a vortex-glass phase occur over a 
significant range x > xvl ? 

In this work, we argue that for rational values of /, 
the superconducting transition temperature of the ar- 
ray vanishes, for increasing dilution, at a critical value 
Xs(f)- The vortex-lattice ordering remains stable up 
to Xvl(I) > %s{f) but both values are much below 
the value x p corresponding to the geometric percolation 
threshold. For xvl < x < x p there is a zero-temperature 
vortex glass. These features are verified numerically for 
/ = 1/2, using a bond-diluted frustrated XY model on 
a triangular lattice, and extensive zero temperature cal- 
culations. Domain- wall energy calculations gives an es- 
timate of a wide range, xvl < x < x p , for a zero- 
temperature vortex glass below the geometrical percola- 
tion threshold x p = 0.652. We find xs = 0.14(1) and Xvl 
= 0.17(1) consistent with the proposed scenario. In the 
vortex-glass phase, Vt ~ 1 -9 , as estimated from the size 
dependence of defect energies excitations. Interestingly 
enough, this estimate is V£ry close to the value obtained 
for the gauge-glass modelEI which may suggest a common 
universality class. 

We consider a two-dimensional Josephson-junction ar- 
ray in a magnetic field B described by the Hamiltonian 
of a frustrated XY model 



9j Ai 



(1) 



where 6i is the phase of the condensate wave function in 
a grain at site i and is the Josephson coupling. The 
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summation is taken over all nearest neighbors of a reg- 
ular reference lattice. The dimensionless line integral of 
the vector potential Aij about each elementary reference- 
lattice plaquette of area S is X) P ^kj = 27r/, where the 
frustration parameter / = BS/& measures the number 
of flux quanta <I> per plaquette. A bond-dilution con- 
centration x corresponds to Jij being zero or J, with 
probabilities x and 1 — x, respectively. Since any closed 
loops of nonzero bonds J have an area which is an integer 
multiple of the elementary area S, the properties of this 
model are periodic in / with period 1, and it is therefore 
sufficient to consider < / < 1. 
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FIG. 1. Schematic phase diagram of a diluted JJA as a 
function of temperature T and dilution x, for an average ra- 
tional frustration /, showing a superconducting phase (S) a 
normal ordered vortex-lattice phase (N) and a vortex-glass 
phase (VG). A single transition is assumed at x = 0. If the 
transitions are already separated at x = 0, dilution would 
further increase the separation. At T = 0, a vanishing linear 
resistance is expected for x < x p due to vortex pinning. 

For / = 0, the Hamiltonian reduces to the standard 
diluted-jXY model, which is known to be superconduct- 
ing foru x < x p . When / =/= 0, there must also be a 
threshold xvl for vortex-lattice disordering!] below the 
percolation threshold, xvl < x p . In the undiluted case 
x = 0, the ground state for rational / = p/q (q > 2) 
consists of a pinned vortex latticed with a q x q unit cell. 
For small dilution x << xvl, the long-range order of 
the vortex lattice persists, provided an infinite cluster of 
these cells exist. Since xvl(J) corresponds roughly to 
the percolation threshold for cells of size qxq, the perco- 
lation threshold for q x q cell dilutation is reached much 
below the unit bond-dilution threshold. Alternatively, 
long-range order of the vortex lattice requires connectiv- 
ity over at least q bonds, as in bootstrap percolationEj, 
which is known to lead to a percolation threshold below 
the unit bond percolation. Since vortex lattipe disorder- 
ing leads to suppression of phase coherenceO, xvl is an 
upper bound for the superconducting threshold x$. This 
implies that the transition temperature should vanish at 
an xs < x p , and that the thresholds are as illustrated in 
Fig. 1. At least for low-order rational values of /, we 
would expect xs(f') < %s(f) if /' < / since /' requires 
a higher connectivity. For xvl(J) < x < x p , there is 



no long-range order, and this phase should correspond 
to a two-dimensional vortex glass, where Ehjttpe phase 
transition is known to occur only at T = 0BO. An in- 
tervening glass phase near percolation threshold is also 
expected from mean field theoryu. This phase can be 
characterized!] by a critical exponent 9 that determines 
how low-energy excitations AE(L) from the ground state 
behave at long length scales L. For a T = vortex 
glass AE ~ L e , with 8 < 0, and thermal excitations 
of scale £ ~ T~ UT destroy the glass order at any finite 
temperature, leading to an identification of the thermal 
correlation length exponent as vt — V|0|> Ovrr numer- 
ical results for / = 1/2, described below, are consistent 
with this behavior, and provide an estimation of vt- In 
absence of thermal fluctuations, at T = 0, vortices are 
pinned by disorder and a nonlinear response to an ap- 
plied current is expected leading to a vanishing linear 
resistance and nonzero critical current for x < x n . 



Ox=0.05 
□ x=0.10 
Ox=0.13 
Ax=0.14 
<x=0.15 
Vx=0.16 
>x=0.17 
+ x=0.18 
xx=0.19 
*x=0.21 
Ox=0.25 




FIG. 2. Finite size behavior of defect energy [Ei] probing 
the superfluid density for increasing dilution x and various 
system sizes L. The change in the L dependence determines 
the threshold xs- 
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FIG. 3. Finite size behavior of defect energy [E2] probing 
the vortex-lattice lattice stability for increasing dilution x and 
various system sizes L. The change in L dependence deter- 
mines the threshold xvl- 
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We have carried out a detailed numerical study for 
/ = 1/2 at zero temperature, using a bond diluted frus- 
trated JJA on a triangular lattice, where the critical 
dilution threshold for bond percolationB is x p — 0.652. 
For this value of /, vortex-lattice ordering can be conve- 
niently described in terms of a Z 2 chirality order param- 
eter x — J2<ij> i®i ~ ®3 — ^u )/(2 7r )i where summation is 
taken about an elementary plaquette of the actual lattice, 
and the gauge-invariant phase difference is restricted to 
the interval [— w, In the undiluted case, the ground 
state consists of a pinned vortex-lattice corresponding 
to an antiferromagnetic arrangement of x = ±1/2. To 
study the stability of the ordered-phases , we use a defect 
energy renormalization analysis^ at T = 0. A defect is 
created in a system of size L x L by imposing a change 
in the boundary conditions in one direction. The change 
AE(L) in the ground state energy for small systems 
is calculated for a large number of samples by directly 
searching for the minimum energy. We used an improved 
algorithm based on Ref. [H| Typically, 3000 configura- 
tions of disorder have been used for each system size. To 
study both phase coherence and vortex-lattice order, we 
consider two types of defects: (i) From the energy differ- 
ence between periodic E p and antiperiodic E a boundary 
conditions in the phases 0i we obtain AE\ = E a — E p , 
which is a measure of phase coherence, and is related to 
the renormalized stiffness constant J(L) — p AEi/2tt 2 , 
where p = 2/V3 is a geometrical factor for the trian- 
gular lattice. In the thermodynamic limit, J is finite in 
the phase coherent state and vanishes in the incoherent 
state; (ii) A domain-wall defect energy is obtained as 
AE 2 = E r — E p , where E r is the .ground state energy 
with reflected boundary conditions^, corresponding to 
the energy cost for a domain wall in the vortex lattice. 
In presence of disorder, AE\ and AE 2 fluctuate between 
samples, with a distribution that can be characterized by 
its moments. Stability of the ground state against ther- 
mal fluctuations requires that the average [AS], where 
j ] denotes a disorder average, is finite or increases with 
L for the U(l) and Z 2 symmetries respectively. Fig. 2 
shows the behavior of the [AI?i] as function of L for 
increasing dilution. For small x , it increases with L, in- 
dicating the existence of long range phase coherenceO. 
For sufficiently large x it clearly decreases for increasing 
L, indicating a disordered phase. The change in the be- 
havior yields an estimate of xs — 0.14(1) . Fig. 3 shows 
a similar plot for [Ai? 2 ] • The increasing trend with L for 
small x corresponds to a vortex-lattice ordered phase, 
which persists for a small but finite range above xs- For 
large x, it decreases with L, and yields an estimate of 
xvl = 0.17(1) . Thus xvl > xs, as indicated in Fig. 1. 
The disordered phase for xvl < x < x p can be regarded 
as a vortex glass, since it lacks long range order in the 
vortex lattice. 

The stability of the glass phase against thermal 
fluctuations is determined by the size dependence of 
the second moment of the energy excitations Wi — 
y/[AEf] - [AE,] 2 oc L e . Here 6 > indicates a glass 
phase at nonzero temperature, whereas 9 < implies 



that arbitrarily low energy excitations at long length 
scales can be thermally excited, destroying the glass 
phase at any finite temperature£3. The size, dependence 
of w, for a value of x = 0.3 in this regiont£l is shown in 
Fig. 4 and clearly indicates a negative 6 for both w\ and 
w 2 , and so the vortex glass only occurs at T = . The 
exponent vt — 1/ | | of the superconducting thermal 
correlation length £ oc T~ VT can be estimated from the 
slope of wi in a loglog plot, giving vt ~ 1.9. Interestingly 
enough, this estimate is very close to the value obtained 
for the gauge-glass modela, suggesting a common univer- 
sality class, but further data would be necessary to check 
whether vt is x- dependent. 
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FIG. 4. Finite size behavior of the second moment of 
the defect energy distribution wi and W2 in the region 
xvl < x < x p for x — 0.3. The negative slope of loguJi 
x logL gives an estimate of 1/vt- 

At finite temperatures, thermally excited vortices and 
disorder effects can significantly reduce the ordered 
phases for x < xvl(J) since bond dilution introduces 
correlated randomness in the flux, as in the case of jan 
array with disorder only in the positions of the grainsEj. 
Unlike positional disorder, random dilution does not ex- 
plicitly affect the phase difference Qi — 8j between two 
superconducting grains in Eq. (1). Its relevance can be 
studied through two coupled frustrated XY models 

H= -fE<«> [cosfA-^-A^) 

+ cos(<^j - <f)j - A i3 - ntij)] 

- A E< ( 2 ) 

where tij is 1 or with probability x and 1 — x respec- 
tively. In the limit h — > oo, the phases are coupled 
6i = 4>i, and the original model in Eq. (1) is recovered. 
The second term has the same form as the Hamiltonian 
describing positional disorder in a superconducting ar- 
ray in the presence of magnetic fielcuj, with a partic- 
ular bimodal distribution of tij. A detailed analysisjp. 
the small h limit combined with known T ^ resultsa-EI 
for x — and the above calculations at T = 0, sug- 
gest the phase diagram of Fig.l. For coupled XY models 
without disordertj, the coupling h renormalizes to large 
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values even when initially small, while the phase transi- 
tions can be described in terms of vortices in the average 
phase variable (Oi + (pi)/ 2. Guided by this, we consider 
initially the two XY models in Eq. (2) to be indepen- 
dent, and consider the particular rational value, / = 1/2, 
where the relevant excitations, chiral .domain walls and 
vortex charges, are better understoodEO. In this case, 
the disorder variables act as random bonds on the chiral 
order parameter x, and as random dipoles on the vortex 
charges. If the transition in the pure case is single (simul- 
taneous disordering of the chiral and XY-like variables), 
the differently-acting disorder can thus separate the two 
transitions, with vortexr-unbinding at temperatures be- 
low the chiral transitionli3. In fact, Monte Carlo simu- 
lations for the frustrated XY model on a square lattice 
with positional disorder, are consistent with the splitting 
into two transitionsEZl. For the triangular lattice consid- 
ered here, we have estimated the chiral transition tem- 
perature at x = xs, where T$ = (Fig. 1), from the 
peak in the chiral susceptibility and found Tvl — 0.27(21 
which can be compared with the estimated separationo 
AT C = 0.01 at x = 0, if one assumes a double transition, 
which clearly shows that disorder tends to separate the 
transitions. The chiral transition is expected to be in 
the universality class of the random bond Ising model, 
where recent studies have shown that the specific heat 
has a broad peak with a very weak log log (T — T c ) di- 
vergence but the other exponents remain with the pure 
Ising model valueaia. This is consistent with Monte Carlo 
simulations of-thc frustrated XY model on a site-diluted 
square latticeO, where it is found that the specific heat 
has a broad peak which does not clearly grow with lat- 
tice size, in contrast to the undiluted case which grows 
almost logarithmically. Even when a finite coupling be- 
tween the two terms in the Hamiltonian of Eq. (2) is 
taken into account, the effects of disorder on the chiral 
order parameter should still remain, since the coupling 
term should essentially lock equivalent vortices and chi- 
ral variables in both phases 9i and <pi ■ For other values of 
/, we expect similar qualitative behavior, as illustrated 
in Fig. 1, but with the chiral transition replaced by the 
thermal disordering transition of a vortex lattice with a 
higher order discrete symmetry. 

Experimentally, the vortex glass phase for xvl < x < 
x p could be identified through the change in the current- 
voltage characteristicsO extracting the critical exponent 
vt- Another signature would be the disappearance of 
ordered-phase resistance- minima at / = p/q when x is 
in the /-insensitive vortex-glass region xvl(I = 1/2) < 
x Xp . 
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